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QUENCHING  OF  ADAPTIVE  CONTROL  SYSTEM  RESPONSE  TO  TEST  SIGNAL 

by 

Rufus  Oldenburger* 
and 

Luther  Schrock** 

Purdue  University 


In  adaptive  control^  test  signal  may  be  used  to  identify  the  parameters  of 
the  system  to  be  controlled.  A-teet  mi gwl  disturbs^ tha_fflraimm.  It  is  therefore 
desirable  to  eliminate  the  effects  of  this  signal  as  soon  as  the  identification 
has  been  complete^^  After  identification  a  quenching  signal  may  be  introduced  to 
eliminate  the  system Vesponse  to  the  test  signal.  The  systems  treated  here  are 
described  by  linear  diA(s£|milH^equationa  with  slowly  varying  coefficients.  The 
nature  of  the  quench  ing'^ifcsmjf'aspends  on  whether  or  not  it  is  bounded  or  unbounded. 
The  unbounded  quenching  signal  is  a  linear  combination  of  a  properly  weighted 
Impulse  and  derivatives  of  an  impulse.  The  weights  are  determined  by  the  initial 
conditions  on  the  system  at  the  instant  of  quenching  and  the  system  parameters. 
Unbounded  quenching  is  considered  to  be  optimum  if  the  response  to  the  test  signal 
is  eliminated  with  minimum  integral  squared  error.  The  bounded  quenching  signal 
is  obtained  by  scheduling  the  lengths  of  time  its  value  is  either  at  the  upper  or 
lower  bound.  The  quenching  signal  is  determined  by  the  test  signal  and  the  system 
to  be  controlled.  Therefore  as  soon  as  the  system  is  identified  quenching  can  be 
accomplished  by  scheduling  regardless  of  the  other  disturbances  to  which  the  system 
is  subject.  The  method  applies  to  the  ouenching  of  system  response  whether  or  not 
adaptive  control  is  involved. 


*  Director,  Automatic  Control  Center;  Professor  of  Mechanical  Engineering 
**  Research  Assistant,  School  of  Mechanical  Engineering 
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INTRODUCTION 

State  space  for  a  system  described  by  a  differential  equation  of  nth  order 
le  defined  as  a  system  of  cartesian  coordinates  where  the  ordinates  are  the 
controlled  variable  and  its  first  n-1  derivatives.  The  normal  state  of  a  system 
is  the  state  due  to  the  normal  input.  The  normal  input  is  t.he  input  to  the  system 
in  the  absence  of  the  test  or  quenching  signal. 

Adaptive  oontrol  is  indicated  where  the  parameters  of  the  system  to  be 
controlled  are  time  varying.  It  is  important  to  know  the  value  of  the  parameters 
to  match  the  controller  to  the  system.  The  determination  of  the  parameters  of 
the  system  is  commonly  referred  to  as  the  "identification  problem".  One  method 
of  identification  is  to  use  the  normal  input  and  the  corresponding  system  response. 
The  other  is  to  introduce  a  test  signal.  The  test  signal  disturbs  the  system. 
The^esnonse  to  the  normal  input  is  often  difficult  to  distinguish  from  the  noise 
present.  In  such  cases  a  test  signal  for  which  the  system  response  dominates  the 
noise  is  preferred  to  the  normal  input. 

Mishkin,\j3raun,  Corbin,  Merriam  and  others  present  technioues  for  solving  the 
identification  ptnjblem  .  1  .  j 

We-prepeee  the  adaptive  control  ehcarfn  irj Figure  1.  The  switches  are  closed 
during  normal  operation  and  open  for  the  period  of  Identification.  The  switches 
may  be  omitted,  but  then  the  computation  becomes  much  more  difficult.  A  test 
signal  is  used  for  identification  to  eliminate  the  effects  of  the  test  signal. 

The  theory  is  restricted  to  controlled  systems  described  by  linear  differential 
equations  with  constant  coefficients.  The  form  of  the  equation  is  known  at  the 
start  but  not  the  values  of  the  coefficients  which  are  parameters.  It  is  assumed 
that  the  parameters  of  the  system  vary  slowly  enough  so  that  they  may  be  considered 
constant  during  identification. ^Immediately  after  identification  the 
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paraaeter  values  are  available  to  calculate  the  ouenching  signal.  The  response 
to  the  test  signal  may  then  be  quenched  by  scheduling  determined  only  by  the 
test  signal  and  system  characteristics.  This  signal  is  independent  of  the  other 
disturbances  to  which  the  system  is  subject. 

The  nature  of  the  quenching  signal  depends  on  whether  or  not  it  is  bounded. 

It  will  be  shown  that  the  unbounded  quenching  signal  is  a  linear  combination  of  a 
properly  weighted  impulse  and  derivatives  of  an  impulse.  The  order  of  the  highest 
derivative  is  one  less  than  the  order  of  the  system.  The  weights  are  determined 
by  the  state  of  the  system  at  the  instant  of  quenching  and  the  system  parameters. 
Unbounded  quenching  is  considered  to  be  optimum  if  the  response  to  the  test  signal 
is  eliminated  with  minimum  integral  squared  error.  The  bounded  quenching  signal 
is  obtained  from  a  "bang-bang**  controller.  The  signal  is  either  positive  or 
negative  and  its  absolute  value  is  less  than  or  equal  to  a  constant.  Oldenburger, 
LaSalle,  Bellman,  Rosonoer  and  others  present  methods  of  obtaining  the  control 
functions  necessary  to  determine  when  the  signal  is  positive  or  negative .2,  3,  U,  5 

The  control  time  is  defined  as  that  time  interval  during  which  the  bounded 
quenching  signal  is  either  positive  or  negative.  The  control  functions  for  a 
bang-bang  controller  are  used  to  obtain  the  control  times.  These  time  intervals 
are  expressed  as  functions  of  the  test  signal,  the  magnitude  of  the  ouenching 
signal  and  the  system  to  be  controlled.  The  quenching  signal  is  introduced  by 
scheduling  the  control  times.  For  some  systems  the  control  times  are  independent 
of  the  system  parameters.  A  computer  may  be  readily  programmed  to  compute  the 
control  times  if  they  are  not  independent  of  the  system  parameters. 

Not  all  systems  are  treated;  however  the  method  is  satisfactory  for  a  large 
class  of  simple  systems.  For  the  general  second  order  and  higher  order  systems 
the  expressions  for  the  oontrol  times  become  complicated  and  impractical  without 
approximations. 
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GHJERAL  THEORY  FOR  CASE  OF  UNBOUNDED  QUENCHING  SIGNAL 


The  general  system  is  shown  in  Figure  2.  It  is  assumed  that  this  system  is 
linear  and  that  g(t)  is  its  impulse  response.  The  Laplace  transform  G(s)  of  g(t) 
represents  the  transfer  function  of  the  system. 

The  signal  r(t)  may  be  the  normal  input  or  an  equivalent  input  due  to 
disturbances  entering  the  system  at  other  points.  The  extra  signal  m(t)  is 
composed  of  two  components  m^(t)  and  ntq(t).  Here  m^(t)  is  the  test  signal  for 
identification  and  «q(t)  is  the  quenching  signal  introduced  to  drive  the  system 
to  the  normal  state.  In  general  r(t)  is  a  function  of  time,  either  deterministic 
or  random. 

Since  the  system  is  assumed  linear,  superposition  holds  and  the  output  c(t) 
la  the  sum  of  the  responses  to  the  three  inputs  r(t),  m^(t)  and  iQq(t).  The  three 
outputs  are  given  by  the  convolution  integrals 


ci(t) 

Cq(t) 

cr(t) 


r « 


I  g(t-<r)mi(r-)d«r 


T)mq(T)dT 
A)r  ( A)dA 


(1) 


where  cr(t),  c^(t)  and  cq(t)  are  the  responses  to  r(t),  m^(t)  and  n^(t)  respec¬ 
tively.  Thus 


c(t)  -  ci(t)  ♦  cq ( t )  ♦  cr(t)  . 


(2) 


We  define  the  system  error  eq  by 


-  c. 


(3) 


where  ct  is  the  actual  output  of  the  system  and  c^  is  the  desired  output. 


When  m^(t)  is  introduced  the  response  to  this  signal  is  superimposed  on  the 
response  to  r(t).  If  the  identification  requires  time  T  and  mi(t)  is  introduced 
at  t  *  0,  it  is  desirable  to  return  the  system  immediately  to  the  normal  state 
at  t  •  T.  Ideally  we  would  like  t>  return  the  system  with  zero  error.  Figure  3 
shows  the  desired  response. 

For  O^taT  the  system  response  is  the  sum  of  cr  and  c^.  For  t*»T  we  then 
want  the  output  cr.  Therefore  we  define  the  desired  output  cd  by 

cd  -  cr  tswT  .  (4) 

We  have  stated  that  ideally  the  system  response  should  be  returned  to  the 
normal  state  with  zero  error;  this  corresponds  to  moving  from  one  point  in  state 
space  to  another  in  zero  time.  Due  to  limitations  on  the  system  and  on  the  power 
available,  a  finite  time  will  be  required  to  return  the  system  to  the  normal  state. 
The  actual  response  of  the  system  is  shown  in  Figure  4> 

Since  the  system  is  not  returned  to  the  normal  state  instantaneously  at 
t  •  T  the  error  e^  is  given  by 

eq  “  ca  "  cd  “  <C1  *  cq  *  cr)  "  cr  "  C1  *  cq  t^T  •  (5) 

Since  Cq(t)  ■  0  for  t^T  it  is  convenient  to  introduce  a  new  variable  t' 

where 

t»  -  t  -  T 

and  replace  Cq(t)  by  a  new  function  c2(t')  where 

c2(t»)  -  Cq(t )  . 
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Equation  (5)  may  be  written  in  terms  of  the  new  variable  t*  as 

e(t')  *  c^(t‘  ♦  T)  +  C2(t')  t'»0  (6) 

where 

e(t')  -  «q(t)  . 


We  introduce  n^t')  where 

m2  ( t 1 )  -  mq(t) 


The  integral  squared  error  I  is  given  by 

I 


-  J  e2(t')  dt'  . 


By  Equations  (6)  and  (7)  we  obtain 

■  OO 


1  "  *  2Cl^'  *  7)c2(t  1 )  *  c|(t ' )]  dt' 

The  error  and  hence  the  integral  squared  error  is  aero  if 


c2(t')  •  -  ci(t'  *  T) 


t'»»0  . 


Now  c2(t')  is  given  by  the  convolution  integral 

r°° 

c2(t')  -  /  g(t'  -T)ra 2(T)  d T  . 

J-OO 

By  Equations  (9)  and  (10) 

r°° 

/  g(t'  -T)m2(T)  d T  -  -  c^t'  ♦  T)  . 


(7) 


(«) 


(9) 


(10) 


(11) 


Equation  (ll)  may  be  solved  for  m2(t')  by  taking  the  Laplace  transform  of 
both  sides  to  obtain 


G(s)M2(s)  -  -  ♦  T)] 


(12) 
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where  H2(s)  is  the  Laplace  transform  of  m2(t')  with  respect  to  the  time  variable 
t  • ,  and  T)j  is  the  Laplace  transform  of  the  output  c^(t'  ♦  T)  with 

respect  to  t1. 

Since  mi(t)  *  0  for  t^T  the  response  of  the  system  does  not  depend  on  m^(t) 
for  t»»T,  but  only  on  the  initial  conditions  at  t  ■  T;  that  is,  t'  *0. 

In  general  the  linear  system  under  study  is  given  by  the  differential  equation 


dn<5l  „  ,  dn‘lcl  «. 

an-l“^T 


dc-i 

♦  _  _  X  ♦ 

al  dt  ao°l 


q  dt<i 


♦  bq.! 


df*-1 


dm 

*  bi 

*  *  1  dt 


“I  ♦ 


b0ml 


(13) 


for  real  coefficients  a0,  a^,  .  .  .  an_i  and  b0,  .  .  .  bq.  We  assume  that  qjj  n-1. 
Por  each  1  we  have 

d^ci  d^ci 


dv 


d(t'  J1 


The  initial  conditions  at  t  »  T  are 


c1(T)  -  cc 


c’(t)  -  c; 


^-l^T)  .  cfr-l) 


where  cQI  c^,  .  .  .  c^  ^  are  the  values  of  c^,  c^,  .  .  .  c^“^  at  t 1  *  0. 

Let  C^(s)  denote  J^£c^(t'  *  T )j  for  the  time  variable  t'.  By  Equation  (13) 

£anC^  —  s^^^Cp  —  s*^~^Cp  —  ...  -  Cq~^  *  an-l[3n~^Ci-sn“^c0-sn_3cp-. . 

*l[sCl  -  co]  4  a0Cl  ■  0 


(H) 
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whence 


We  note  that  G(e)  has  the  aane  denominator  as  C^(s)  but  that  the  numerator 
of  G(a)  la 

V  •  Vi*""1  -  ■  ■  •  *  »0  • 

By  Equation  (12) 


Now  C2(t • )  ia  the  aolution  of  the  differential  equation  witn  initial 
eonditiona 

C2(0)  -  -  ci(T) 

c2(0)  •  -  c^(T) 


(^-1^0)  -  - 

Since  c2  ■  0  for  t  0  the  problem  ia  to  create  the  above  initial  valuee  of 
c2  and  ita  derivatives. 

If  the  inverse  Laplace  transform  of  both  sides  of  Equation  (16)  ia  taken, 
then  m2(t')  will  be  that  signal  which  produces  the  negative  of  the  response  c^(t) 
for  t»T.  Figure  5  shows  the  response  c2(t'). 

If  c{(T)  /  0  the  response  C2(t’)  requires  an  instantaneous  change  in  velocity 
and  hence  an  infinite  acceleration  of  the  system  output.  For  a  physical  system 
there  is  always  mass  or  inertia  associated  with  the  output;  therefore  Infinite 
acceleration  reauires  infinite  force. 


I 
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SECCND  AND  THIRD  CRDER  CASES  FOR  UNBOUNDED  QUENCHING  SIGNAL 


We  consider  a  second  order  system  with  the  transfer  function 

kl 


G(s) 


♦ 


ais  ♦  ae 


(17) 


The  differential  equation  describing  the  system  la 


d^c-i  dci 

- -  „  ♦  ai  — —  *  a0ci  -  0 

d(t')2  dt« 


where  the  input  is  aero  since  we  consider  only  the  response  for  t'^O  and 
m^(t)  •  0  for  t' »  0. 

Taking  the  Laplace  transform  of  both  sides  of  Equation  (18)  yields 
(s2  ♦  aja  ♦  ae)C1(s)  -  (ac0  ♦  c£  ♦  a^)  ■  0 

where  c0  and  c,J  are  the  initial  conditions  at  t*  *0. 

Solving  Equation  (19)  for  C^s)  we  obtain 

^(s)  -  e0B  ♦  alCp  ♦  Cp 


(18) 


(19) 


s2  ♦  a-^s  ♦  a0 

Substituting  C^(s)  from  Eouation  (20)  and  G(s)  from  Equation  (17)  into 
Equation  (12)  and  solving  for  1^(8 )  we  get 


(20) 


M2(s)  "  "  “k^  [c°“  4  (aic°  4  C<^] 


(21) 


The  unit  impulse  6(t')  may  be  defined  by 


S(  f) 


lira 

a^O 


u(t' )  -  u(t'  -  a) 


(22) 
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where  u(t')  and  u(t*  -  a)  are  unit  step  Amotions  at  t '  "0  and  t'  ■  a 
respectively.  The  derivative  of  a  unit  impulse  is  taken  as 


d  S  (tQ 

dt' 


lim 

a-^0 


u(t')  -  2u(t'  -  a)  ♦  u(t'  -  2a) 

- 77 - 


(23) 


in  the  sense  that  the  limit  of  the  Laplace  transform  of  the  fraction  in 
Equation  (23)  as  a^O  is  s.  The  plot  of  this  fraction  versus  t'  is  a  double 
pulse,  three  of  vhich  are  pictured  in  Figure  7c. 

Similarly  the  second  derivative  is  taken  as 

d2l(t')  .  lim  u(t 1 )  -  3u(t'  -  a)  ♦  3u(t 1  -  2a)  -  u(t’  -  3a)  (2a) 

dt'  a-*0  a3 

where  the  limit  of  the  Laplace  transform  of  the  fraction  in  Eouation  (24)  is  a2. 
The  plot  of  this  fraction  is  a  triple  pulse. 

We  shall  illustrate  the  use  of  the  fraction  on  the  right  in  relation  (22). 
Consider  a  simple  second  order  system  described  by  the  differential  equation 

c"  -  kx  m2  (t * )  (25) 

where  the  primes  denote  differentiation  with  respect  to  t ' . 

Let  C(s)  be  the  Laplace  transform  of  c(t')  with  respect  to  the  time  variable 
t*  and  M2(s)  be  the  Laplace  transform  of  n^t')-  Taking  the  Laplace  transform  of 
both  sides  of  Equation  (25)  and  solving  for  C(s)  we  obtain 

c(s )  ■  k-iMsia)  •  se° ;  .g° 

where  c  ■  c0  and  c'  ■  Cq  at  t 1  *  0. 


(26) 


-li¬ 


the  approximation  to  a  unit  impulse  and  its  derivative  is  given  by 
Equations  (22)  and  (23)  when  a  /  0.  We  let 

a2(f)  ■  -wx  -au-(tl  -  al]  -  w2  [  1 ).  ■-  Mi  P.L*'  _.2a)]  (27) 

where  is  the  weight  of  the  approximate  impulse  and  W2  is  the  weight  of  the 
approximate  derivative. 

Taking  the  Laplace  transform  of  both  sides  of  Equation  (27)  we  get 

.-as  -2aa 


«.<■)  •  -<?•»*  * 


(28) 


Substituting  M2(s)  from  Equation  (28)  into  Equation  (26),  taking  the  inverse 
Laplace  transform  and  collecting  t  ms  we  obtain 


c(t') 

c'(t') 


kjW^a  _  k^Wj  —  kjWjt 1  ♦  c0  ♦  c^t 1 


c©  -  kxWx 


(29) 


We  now  let  a-eO  in  Equations  (29).  The  result  is 
c(t')  •  ce  ♦  Opt*  -  k1(W2  ♦  Wxt • ) 


c'(t') 


co  *  klwi 


(30) 


If  c(t')  •  c'(t')  ■  0  for  t'*0  we  must  have 


W, 


Wo 


*1 


(31) 


Letting  a-eO  in  Equations  (29)  gives  the  same  result  as  letting  a-*0  in 
Equations  (22)  and  (23). 
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At  t'  ■  0"  we  have  c  -  cQ  and  c'  -  Cq.  At  t'  -  0*  we  have  c  -  c*  -  0. 
Therefore  ce  and  are  reduced  to  zero  Instantaneously  at  t'  *  0  for  the  values 

of  and  W2  given  by  Equations  ( 31 ) . 

We  now  take  the  inverse  Laplace  transform  of  both  sides  of  Equation  (21) 
to  obtain 

"  "  "kj  [co<tf8>  4  <alco  4  OKf)]  (32) 


»diere  €( t')  is  the  unit  impulse  and  ^(s)  is  the  derivative  of  the  unit  impulse, 
referred  to  as  the  "doublet'.'  6  ,  7.  The  initial  condition  c0  is  the  weight 
associated  with  the  doublet  and  (a^c0  ♦  Cq)  is  that  of  the  impulse. 

If  the  system  is  third  order  with  the  transfer  function 


0(a) 


kl 

2 


.3 


•2*  ♦  »!» 


(33) 


M2(s)  is  given  by 


M2(s) 


(®o  4  •2°o)>  4  (co  4  *2co  4  *100)]  . 


(34) 


The  inverse  Laplace  transform  of  both  sides  of  Equation  (34)  now  yields 

i»2(t')  •  -  ^  [c0£(e2)  ♦  (cp  ♦  a2c0)X(a)  ♦.  (cj  ♦  a2Co  ♦  •ic0)4(t')]  (35) 

where  /V)  is  the  second  derivative  of  the  unit  impulse  and  is  referred  to  as 
the  "triplet".  The  weights  of  the  triplet,  doublet  and  impulse  are  cG,  (c^*a2c0) 
and  (c^ajc^ajc^  respectively. 

It  Is  noted  that  the  weights  of  the  impulse  and  its  derivatives  are  deter¬ 
mined  from  the  initial  conditions  at  t'  “0  and  the  system  parameters.  The  order 
of  the  highest  derivative  is  one  less  than  the  order  of  the  system. 
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The  sigial  n^t')  ia  not  physically  realizable  due  to  the  appearance  of  the 
Impulse  and  its  derivatives  whence  the  error  e(t*)  will  not  be  zero. 

Truxal  shows  that  the  optimum  manner  for  an  ideal  system  to  reach  the  origin 
from  some  point  in  the  phase  plane  is  as  shown  in  Figure  6.  8.  The  path  begins 
at  some  initial  condition  c0,  Cp  and  goes  to  infinity  then  returns  along  the 
c'  ordinate  from  infinity  to  the  origin.  Since  the  area  under  the  reciprocal  plot 
of  1/c'  versus  c  is  zero,  and  time  in  the  phase  plane  is  given  by 


no  time  is  required  to  reduce  c  to  zero. 

The  unit  impulse  may  be  approximated  by  considering  a  pulse  of  finite  height 
and  short  time  duration  such  that  the  area  under  the  pulse  is  unity. 

Figure  7  a,b  and  c  shows  the  family  of  paths  in  the  phase  plane,  time  domain 
and  the  respective  quenching  signals  for  a  second  order  system.  As  the  pulses 
become  hi^ier  and  of  shorter  time  duration,  the  time  to  reach  the  origin  from  the 
initial  condition  at  t  ■  T  decreases. 

To  introduce  the  properly  weighted  pulses  the  computer  must  determine  c0, 

Cq,  .  .  .at  the  instant  t  ■  T.  In  practice  we  cannot  compute  the  third  or  higher 
order  derivatives  because  of  noise  present  in  the  system.  We  can  however  calculate 
the  initial  conditions  from  the  known  test  signal  and  the  system  equation.  Thus 
the  control  time  may  be  computed  in  terms  of  the  pulse  height  and  the  initial 


conditions  at  t  ■  T. 


SYSTi^M  WITH  BOUNDED  QUU'iCHING  SIGNAL 


The  differential  equation  describing  the  linear  system  under  study  will  be 
taken  to  be 


dnc 

d(t')n 


dn“^c 

n-x  d(f  r-1 


♦ 


dc 

dt' 


aoc 


<lm2 


(36) 


where  Im^  2k?.  Then  the  quenching  signal  is  bounded.  Oldenburger  has  shown 
that  under  rather  general  conditions  the  best  return  to  "equilibrium"  is  attained 
by  operating  "bang-bang" .2 .  For  the  problem  at  hand  equilibrium  is  the  normal 
state  of  the  system.  In  every  reasonable  sense  the  trajectory  is  optimum,  i.e. 
least  time  to  equilibrium,  minimum  overswing  or  underswing  and  minimum  area 
between  the  trajectory  and  the  t*  axis,  etc. 

We  shall  now  make  use  of  the  control  functions  for  a  bang-bang  controller  to 
obtain  the  control  times.  The  quenching  signal,  scheduled  according  to  the  control 
times,  will  bring  the  system  to  equilibrium. 

Example  1;  Slmnle  Second  Order  System 

'We  consider  the  system  described  by  the  differential  equation 

c"  «  kjm2  (37) 

where  the  primes  denote  differentiation  with  respect  to  t'.  The  problem  is  to 
determine  the  control  times  in  terms  of  the  initial  conditions  at  t1  *0. 

The  optimum  control  function  C  for  the  system  of  Lauation  (37)  Is  given 


(3*> 


by  .2 


>c't  c1 
2k}k2 
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The  schedule  for  optimum  control  is 


n»2  *  “  (sgh  E)  k2  (39) 

where  sgn  £  denotes  the  sign  of  £.  Suppose  that  E  /  0  at  t'  *  0,  The  transient 
from  t'  ■  0  to  the  instant  t^  where  E  *  0  is  the  first  phase  of  the  solution. 

This  is  followed  by  a  second  phase  terminating  in  equilibrium.  The  duration  of 

i 

this  phase  will  be  denoted  by  t2. 

By  Equation  (37) 

i 

c'  ■  cD  *  kjH^t' 

C  ■  co  4  cot'  *  (t')2 

where  c  *  cQ  and  c'  ■  c^  at  t'  *  0. 

We  substitute  c'  and  c  from  Equations  (40)  into  £  ■  0  to  obtain 


co  *  * 


:i)2  •  <•*  •■)  (c°  ■  ° 


where  sgn  c'  denotes  the  sign  of  c'  when  E  ■  0. 

Solving  for  t{  from  Equation  (41)  we  obtain 

t‘  -  (sgnE)cp  ♦V|(cp)2  ♦  (sgnE)kik2CQ  ^  ^ 

1  *1  k2 

Let  c^  and  be  the  values  of  c  and  c'  at  t'  "  tj_.  By  Equat  ions  (40)  and 


(42)  we  have 


c{  ■  -  (sgn  E)  V  i(c©)2  ♦  (sgnt)k1k2Cp 


The  initial  conditions  for  the  last  phase  are  and  c^.  Letting  t"  ■  t'-t^ 


for  the  last  phase  we  have 


♦  kjn^t" 
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I4ien  t"  •  the  system  la  at  equilibrium  and  c'  ■  0.  Thus 


H 


-  C1 
kl®2 


(45) 


Substitution  of  from  Equation  (43)  into  Equation  (45)  yields 

t‘  -  V  *(co)Z  4  (a*n  C  )klk2co  .  (46) 

kl  k2 

Equations  (42)  and  (46)  give  t^  and  tg  respectively  in  terms  of  the  initial 
conditions  c0  and  c^.  It  is  unlikely  that  the  case  where  I  ■  0  »t  t  ■  T  will 
occur.  This  caee  is  taken  care  of  by  having  t£  ■  0.  The  times  tj  and  t£  approach 
sero  as  kj  becomes  arbitrarily  large. 

Equation  (37)  is  the  differential  equation  of  the  syetem  after  identification. 
We  let  »i(t)  be  the  test  signal,  to  be  introduced  at  t  ■  0.  The  differential 
equation  of  the  system  for  t'««0  is 


c" 


kl  "1 


(47) 


where  the  primes  denote  differentiation  with  respect  to  t. 

We  wish  to  express  t{  and  tg  in  terms  of  the  known  test  signal  and  the 
identification  period  T.  We  let  the  test  signal  be  an  impulse  of  weight  W.  Here 
is  not  bounded  but  it  is  understood  that  the  quenching  signal  is.  The  solution 
to  Equation  (47)  is 


c(t )  ■  kxWt 

c  * ( t )  -  kjW  . 


(4«») 


At  the  end  of  the  time  interval  T  the  Identification  is  assumed  to  be 
complete  and  Equation  (37)  applies.  Thus  the  initial  conditions  for  Equation 
(37)  are 


-17- 


c(T)  -  c0  -  kxWT 

c'(T)  -  Cq  -  kjW  . 


(49) 


Here  cn  and  c' 
o  0 

Into  Equations  (42) 
are  giv.  by 


are  both  positive.  Hence  0.  We  substitute  c0  and 

and  (46)  to  obtain  t|  and  t2  respectively.  The  control  times 


t 

tl 


W  ♦  k2WT 

k2 


(50) 


t 


2 


V  ♦  k2wr' 

k2 


It  is  of  interest  to  note  that 

^2  "  ^1  “  ^/k2  • 

Physically  one  cannot  obtain  an  impulse  but  must  approximate  the  impulse 
by  a  pulse  of  finite  height  and  short  time  duration.  Figure  6  shows  an  approxi¬ 
mation  to  an  impulse  where  H  is  the  height  and  a  is  the  duration.  It  is  assumed 
that  a«*T. 

If  ra^  is  the  pulse  shown  in  Figure  8,  the  initial  conditions  for  Equation 
(37)  are  obtained  by  solving  Equation  (47).  These  initial  conditions  are 

c0  "  fkjH(2aT  -  a?) 

,  (51) 

cQ  -  kjHa  . 

Substitution  of  the  initial  conditions  from  Equations  (51)  into  Equations 


(42)  and  (46)  yields 
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,  Ha  ♦  V$H2a2  ♦  ik^UaT  -  a2)' 
li  '  k2 

,  V} H2a2  ♦  fkpH(2aT  -  a2)' 

‘2  k2 


(52) 


We  define 

H  -  W/a  • 

Letting  a-*0  while  W  remains  constant  Equations  (52)  reduce  to  Eouations  (50). 

If  is  a  step  input  of  hei#it  H  and  duration  T,  the  control  times  become 

t  «  _  HT  ♦  V  ^T2  *  ikgHT2  ' 

(53) 


k2 


t'  -  *  ik2HT2' 


We  now  have  expressions  for  the  control  times  in  terms  of  the  test  signal 
m^,  T  and  kj.  We  note  that  the  control  times  for  the  simple  second  order  system 
are  Independent  of  the  system  and  are  fixed  by  the  test  signal.  Thus  it  is 
relatively  easy  to  schedule  the  quenching  signal  when  the  test  signal  is  known. 


Example  2:  Third  Order  System 

We  consider  the  system  described  by  the  differential  equation 

Tc'"  ♦  c"  -  kjmj  (50 

where  In^l  ^  kj,T  is  the  time  constant,  of  the  system  and  the  primes  denote 
differentiation  with  respect  to  t'. 
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Por  the  system  of  Eauation  (54)  Oldenburger  Introduces  the  two  control 
functions  £  i  and  £  2  where 


Ts  1  -  V  4  -  (8^)klk^ln^lV^l+(sgr^0~^]exp(^!)  } 


,2  -V 


«,  1^ 

2kik2 


(55) 

(56) 


where 

-  c  ♦  T  c  .  (57) 

For  disturbances  normally  encountered  the  log  term  of  £  ^  is  small  compared 
to  the  rest  of  the  terms,  and  can  therefore  be  dropped.  2.  We  therefore  take 

E  ■  E2. 

Optimum  control  is  obtained  by  using  Equation  (39)  end  normally  involves 
three  phases,  for  each  of  which  m2  is  k2  or  -k2» 

We  begin  with  a  set  of  initial  conditions  c0,  c©  and  c©  and  let  the  system 
travel  over  the  first  phase  to  the  first  switch  point  where  £  ■  0.  IXiring  this 
phase  Equation  (54)  becomes 

Tc"'  ♦  c"  ■  -  (sgn  E  )kjk2  .  (58) 

The  solution  to  Equation  (58)  is 

c  -  7*  2£c©  ♦  (sgn£)k2k2j(e-Tt'  -l)  ♦  c©  *  (c©  ♦TCo)t' 

♦  (sgn£  )Tkjk2t '  -  (sgn  E  )fkjk2(t '  )2 


c»  -  (c©  4Tco)  4  (sgn  E  )T  kik2  -t[co  4  (sgn  I  )k]k2]  e“  ^rt ' 

-  (sgn  E)kjk2t  ■ 

c"  -  ^c©  ♦  (sgn  E)k]k2  j  e”Pt'  -  (sgn£)kik2  (59) 
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where  c  ■  c0,  c*  *  and  c"  *  c^  at  t'  ■  0.  We  let  c  ■  c^,  c'  *  and  cM  • 
at  t'  ■  t{;  these  values  of  c,  e*  and  c"  are  the  initial  conditions  for  the  next 
phase.  At  t'  ■  t^  we  form  the  function  y/  where 


W  "  C1 4  r«; 


(60) 


Substitution  of  c^,  c^  and  c^  from  Equations  (59)  into  Equation  (60)  and  the 
equation  for  \f/'  yields 

V',  -  (t{  4  T  )  -  (sgn  r)ik1k2(t{)2  ♦  (c0  -  T2Cq) 

,  ,  (6l) 

¥{'-  U/‘-  (sgnDk^ti 


where  y/  and  y/t  are  the  values  of  yf  and  yj*  at  t*  “  0. 

At  t'  ■  t|,  yt  and  ^/'satisfy  E  ■  0.  Hence 

][]“  Vo*  *T  ^  “  (sgnCJ^kjft^)2  ♦  (c0  -  r  2°o) 

.  .  0 

1, 


2kjk2 


Solving  for  t|  from  Equation  (62)  we  get 


t,  .  (wr)y/'  (sgnpk^1 

1  klk2 


(62) 


(63) 


For  the  second  phase  we  have 

7*  c"'  ♦  c"  •  ♦  (sgn  E)kjk2  . 

The  initial  conditions  are  ci,  c-[  and  c£  given  by  Equations  (59)  at  t'  ■  t^. 
Letting  t"  •  t'-t|  the  solution  of  Equation  (64)  is 


(64) 
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c  -  T2[c£  "  (sgn  E)k]k2  J  (e”^1"  -l)  +  cx  ♦  (c{  ♦  T  ciH" 

-  (sgn  C )t ^1^2!"  4  (sgnZ  )^kik2(t")2 

c'  -  (ci[  ,fTci')  -  (sgn  EJrkjkj  -t[c1  “  (sgn.Cjk^j  e-^1" 

♦  (sgn  C)k^k2t"  (65) 


c"  »  -  (sgnEjk^j  e“7t'"  +  (sgnTjk^  . 


We  let  c  ■  C2»  c'  *  c^,  c"  m  c£  at  t"  ■  t^.  Let and  l^/'  be  the  values  of 
Iff  and  If/'at  t"  •  t^.  Substituting  C2,  and  c!J  from  Equations  (65)  into  If/* 
and  Vi* «  and  substituting  the  result  into  E  “  0,  we  obtain 


yf'(t2  *  T  )  +  (sgn  E)|k1k2(t2)2  ♦  (cx  -  T  2c[) 

♦  (sgnli/*)  (agnpkjkgt^]2 

*  2klk2 


(66) 


Solving  for  tp  from  Equation  (66)  we  obtain 

,  -(sgn  E)W  +  V»(VP2  -  (sgnDki^Vf 


kxk2 


(67) 


We  may  determine  \JJ  and  \J/'  in  terms  of  cD,  c^,  and  t^  by  substituting 
ci,  cj  and  c^  from  Equations  (59)  into  Equation  (60)  and  the  equation  for  \jj*  . 

For  the  last  phase  let  t"'  ■  t1  -  (t^+t;?).  For  this  phase  we  have  Equation 
(5®)  valid  where  the  initial  conditions  are  now  02,  C2  and  Cp,  obtained  from 
Equations  (65)  at  t"  *  tp.  The  form  of  the  solution  for  this  phase  is  given  by 
Equations  (59).  At  t"'  ■  t^  we  have  Cj  ■  cj  ■  c'j  ■  0.  We  may  determine  tj  by 
setting  any  of  the  three  quantities  equal  to  zero.  It  is  sufficient  to  write 


c3  "  [c2  4  (8Cn  C)kik2  j  e”^t3  -  (sgnE)kik2 


(68) 


-  0 
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whence  „ 

t’  -Tln[l.  UnEij^]  .  (69) 

We  may  determine  in  terms  of  cQ,  c^,  c^,  t^  and  ^  by  substituting  c^  from  the 
last  of  Equations  (59)  into  the  last  of  Equations  (65). 

Equation  (54)  is  the  differential  equation  of  the  system  after  identification. 
We  let  m^(t)  be  the  test  signal,  to  be  introduced  at  t  •  0.  The  differential 
equation  of  the  system  for  t'^0  is 

jc"'  ♦  c"  -  kjmx  (70) 


where  the  primes  denote  differentiation  with  respect  to  t. 

Let  the  test  signal  be  an  impulse  of  weight  W.  The  quenching  signal  m2  is 
still  assumed  bounded.  The  solution  to  Equation  (?0)  is 


c(t)  -  k1W(Te“Vt  ♦  t  -T) 


c'(t)  -  kjW(l  -  e-^) 


c"(t) 


klW  -it 


(71) 


At  the  end  of  the  time  interval  T  the  identification  is  assumed  to  be 
complete  and  Equation  (54)  applies.  The  initial  conditions  for  Eouation  (54) 

are 

c(T)  -  cr 


-  k^Te-t7  ♦  T  -T ) 


c'(T)  -  c^  -  kxW(l  -  e-*T) 


'(T) 


-  klW  -  »T 

T'6  * 


(72) 
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Frora  Equations  (72),  y0  and  l^/'0  we  see  that  I>0  at  t'  *  0.  Substitution 


of  c0,  C;  and  from  Eouations  (72)  into  Eouation  (63)  yields 


tl  ’ 


W  ♦  V  iw2  ♦  kpWT 


(73) 


We  now  substitute  cc,  Cq  and  from  Equations  (72)  into  Equations  (59)  to 
obtain  c^,  and  and  substitute  the  result  into  Equation  (60)  to  obtain y/ ^ 
and  if/ [•  Substituting  the  resulting  expressions  for  l^/ ^  and  Into 
Equation  (67)  we  get 


[(k^-W)  ♦  -^i( W-k2t^)2  ♦  k2(|k2(t{)2  -  W(t^T))' j  . 


(74) 


Substituting  cj  from  the  last  of  Equations  (59)  *nto  the  last  of  Eouations 
(65)  and  substituting  the  result  into  Equation  (69)  we  obtain 


t'  -  T  ln<  2  ♦ 


In |2  ♦  -L-  [(£  e"  +  T  ♦  kjje"*^  -  2k2]  e”*1^  j 


(75) 


If  the  impulse  is  approximated  by  a  pulse  of  height  H  and  duration  a,  the 
control  times  become 

,  Ha  ♦  V$H2a2  ♦  kJt*(T  -  fa)' 

h  -  _ 

t2  "  -*2t\2  *  "^^(Ha-k2t^)2  ♦  k2  (*k2(t-[)2  -  Ha(t^T-*a)|'  (76) 

t^  -  T  In ^  2  ♦  -i-  phe~^T(e*#  -l)*k2)  e~^tl  _2k2  e"*t2  j  . 


If  the  teat  signal  ie  a  step  of  height  H  and  duration  T,  the  control  times 


are 
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HT  ♦  V Jh2T2  ♦  ^kjHT2' 


t^  -  T  In  1 2  ♦  ^(H(l-e“^T)*k2)  e"  ▼  L1  -2k2  e~?u-  J  . 


(77) 


We  note  that  t^  and  t2  are  Independent  of  the  system  constants  k^  and  T 
for  all  three  test  signals.  The  first  two  control  tines  are  then  constant  and 
only  t^  muet  be  computed  af ter  T  is  identified.  The  value  of  tj^  and  t2  can  be 
stored  in  the  computer  memory  and  used  to  obtain  t^  immediately  after  identifi¬ 
cation  . 

Higher  order  systems  were  not  considered,  as  the  expressions  for  the  control 
times  become  excessively  complicated.  Also,  if  the  output  is  measured  to 
determine  the  initial  conditions,  derivatives  of  order  hif^ier  than  the  second 
are  difficult  to  obtain  because  of  noise  present  in  the  system. 

For  many  systems  one  cannot  solve  for  the  control  times  explicitly.  This 
ia  true  for  a  general  second  order  system. 

Experimental  Results 

The  results  obtained  for  Example  1  were  verified  on  the  analog  computer. 
Figure  9  shows  the  actual  and  theoretical  system  response  to  the  bounded  quenching 
signal.  The  initial  conditions  on  the  Bystem  at  t’  *  0  correspond  to  the 
conditions  after  identification.  The  control  times  were  calculated  from 
Equations  (42)  and  (46).  It  is  seen  that  the  actual  response  ie  very  close  to 
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the  theoretical.  A  relay  was  used  to  obtain  the  quenching  signal.  The  relay 
had  a  small  amount  of  deadband,  which  accounts  ftr  some  of  the  deviation  from 
the  theoretical  response. 

We  have  thus  found  that  it  is  possible  to  determine  the  control  times  and 
obtain  proper  switching  with  sufficient  accuracy  to  obtain  a  trajectory  which 
is  nearly  the  same  as  the  theoretical  trajectory. 
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Washington  25,  D.  C.  (l) 

Director,  Naval  Research  Laboratory 
Technical  Information  Officer/Code  2000/ 
Washington  25,  D.  C.  (6) 

Commanding  Officer,  Office  of  Naval  Research 

Navy  #100,  Fleet  Post  Office 

New  York,  New  York  (10) 

Commanding  Officer,  0  N  R  Branch  Office 
346  Broadway 

New  York  13,  New  York  (l) 

Commanding  Officer,  0  N  R  Branch  Office 
495  Summer  Street 

Boston  10,  Massachusetts  (l) 

Office  of  Technical  Services 
Technical  Reports  Section 
Department  of  Commerce 

Washington  25,  D.  C.  (l) 

Bureau  of  Ships 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  Code  671  NTDS  (1) 

Bureau  of  Naval  Weapons 
Department  of  the  Nnvv 
Washington  25,  D.  C. 

Attn:  RAAV  Avionics  Division  (l) 


Bureau  of  Naval  Weapons 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  RMWC  Missile  Weapons  Control  Div.  (l) 

Bureau  of  Naval  Weapons 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  RUDC  ASW  Detection  &  Control  Div.  (l) 


Bureau  of  Ships 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  Communications  Branch  Code  686  (l) 

Naval  Ordnance  Laboratory 
White  Oaks 

Silver  Spring  19,  Maryland 

Attn:  Technical  Library  (l) 

David  Taylor  Model  Basin 
Washington  7,  D.  C. 

Attn:  Technical  Library  (l) 

Naval  Electronics  Laboratory 
San  Diego  52,  California 

Attn:  Technical  Library  (l) 

University  of  Illinois 
Control  Systems  Laboratory 
Urbana ,  Illinois 

Attn:  D.  Alpert  (l) 


Air  Force  Cambridge  Research  Center 
Laurence  C.  Hanscom  Field 
Bedford,  Massachusetts 
Attn:  Electronic  Research  Directorate 


Library  (l) 

Technical  Information  Officer 
US  Army  Signal  Research  A  Dev.  Lab. 

Fort  Monmouth,  New  Jersey 

Attn:  Data  Equipment  Branch  (1) 

National  Security  Agency 

Fort  Geo.  G.  Meade,  Maryland 

Attn:  Howard  Campaigne  (3) 

U.  S.  Naval  Weapons  Laboratory  ^ 

Dahlgren,  Virginia  (1/ 


Attn:  Head,  Compution  Div.,  G.H.  Gleissner 
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National  Bureau  of  Standards 
Washington  25,  D.  C. 

Attn:  Dr.  S.  N.  Alexander  (l) 

Aberdeen  Proving  Ground,  BRL 

Aberdeen  Proving  Ground,  Maryland 

Attn:  J.  H.  Giese,  Chief  Compution  Lab.  (l) 

Office  of  Naval  Research 
Resident  Representative 
Purdue  University 
223  Executive  Bldg. 

West  Lafayette,  Indiana  (l) 

Commanding  Officer 
0  N  R,  Branch  Office 
John  Crerar  Library  Bldg. 

86  East  Randolph  Street 


Chicago  1,  Illinois  (l) 

Commanding  Officer 
0  N  R,  Branch  Office 
1030  E  Green  Street 

Pasadena,  California  (l) 

Commanding  Officer 
0  N  R,  Branch  Office 
1000  Greary  Street 

San  Francisco  9,  California  (l) 

National  Bureau  of  Standards 
Washington  25,  D.  C. 

Attn:  Mr.  R.  D.  Elbourn  (1) 

Dynamic  Analysis  and  Control  Laboratory 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts 

Attn:  D.  W.  Baumann  (l) 


Syracuse  University 

Electrical  Engineering  Department 

Syracuse  10,  New  York 

Attn:  Dr.  Stanford  Goldman  (l) 

Communications  Sciences  Laboratory 
University  of  Michigan 
180  Frieze  Building 
Ann  Arbor,  Michigan 

Attn:  Gordon  E.  Peterson  (l) 

Stanford  University 
Stanford,  California 
Attn:  Electronics  Lab., 

Professor  John  G.  Linvill  (l) 


Stanford  University 
Stanford,  California 
Attn:  Electronics  Lab., 

Professor  Gene  Franklin  (1) 

University  of  Illinois 
Urbana,  Illinois 

Attn:  Electrical  Engineering  Dept. 

Professor  H.  Von  Foerster  (l) 

University  of  California 
Institute  of  Engineering  Research 
Berkeley  4,  California 

Attn:  Professor  A.  J.  Thomasian  (1 ] 

University  of  California  -  LA 
Los  Angeles  24,  California 
Attn:  Dept,  of  Engineering, 

Professor  Gerald  Estrin  (1] 

Naval  Research  Laboratory 
Washington  25,  D.  C. 

Attn:  Security  Systems 

Code  5266,  Mr.  G.  Abraham  (1 

Zator  Company 

140  1/2  Mt.  Auburn 

Cambridge  38,  Massachusetts 

Attn:  R.  J.  Solomonoff  (1 

NASA 

Goddard  Space  Flight  Center 
Washington  2?,  D.  C. 

Attn:  Arthur  Shapiro  (l 

Dr.  A.  M.  Uttley 

National  Physical  Laboratory 

Teddington,  Middlesex 

England  (l 

Diamond  Ordnance  Fuze  Laboratory 
Connecticut  Ave.  &  Van  Ness  St. 

Washington  25,  B.  C. 

0RDTL-012 ,  E.  W.  Channel  (l 

Harvard  University 

Cambridge,  Massachusetts 

Attn:  School  of  Applied  Science 

Dean  Harvey  Brook  (l 

Wright  Air  Development  Division 
Electronic  Technology  Laboratory 
Wright  Patterson  Air  Force  Base,  Ohio 
Attn:  WWRNEB  Computer  &  Bionics  Branch  (] 


Laboratory  for  Electronics,  Inc. 

1079  Commonwealth  Ave. 

Boston  15,  Massachusetts 

Attn*  Dr.  H.  Puller  (l) 

Stanford  Research  Institute 
Computer  Laboratory 
Menio  Park,  California 

Attn:  H,  D.  Crane  (l) 

General  Electric  Co. 

Schenectady  5,  N.  Y. 

Attn:  Library,  L.  M.  E.  Dept., 

Bldg.  28-501  (1) 

The  Rand  Corp. 

1700  Main  Street 

Santa  Monica,  California 

Attn:  Numerical  Analysis  Dept. 

Willis  H.  Ware  (l) 

Carnegie  Institute  of  Technology 

Pittsburgh,  Pennsylvania 

Attn:  Director,  Computation  Center 

Alan  J.  Perils  (l) 

Rome  Air  Development  Center,  RCOR 

DCS/ Operations ,  USAF 

Griffiss  Air  Force  Base,  New  York 

Attn:  Irving  J.  Gabelman  (l) 

Air  Force  Office  of  Scientific  Research 
Washington  25,  D.  C. 

Attn:  Dr.  Harold  Wooster  (l) 

Stanford  Research  Institute 
Menlo  Park,  California 
Attn:  Dr.  Charles  Rosen 

Applied  Physics  Group  (l) 

Aeronautics  Division 
Ford  Motor  Company 

Ford  Road,  Newport  Beach,  California 

Attn:  C.  L.  Wanlass  (l) 

Air  Force  Research  Division,  ARDC 
Computer  A  Mathematical  Science  Lab.,  ARCRC 
L.  G.  Hanscom  Field,  CRBB,  Bedford,  Mass. 
Attn:  Dr.  H.  H.  Zschirni  (l) 

Office  of  Chief  Signal  Officer 
Department  of  the  Army 
Washington,  D.  C, 

Attn:  RAD  Division  SIGRO-6D 
Mr.  L.  H.  Geiger 


Bell  Telephone  Laboratories 
Murray  Hill  Laboratory 
Murray  Hill,  New  Jersey 
Attn:  Dr.  Edward  F.  Moore 

Carnegie  Institute  of  Technology 
Graduate  School  of  Industrial  Admin. 
Pittsburgh  13,  Pennsylvania 
Attn:  Dr.  Allen  Newell 

Department  of  Commerce 
U.  S.  Patent  Office 
Washington  25,  D.  C. 

Attn:  Mr.  Herbert  R.  Koller 

Purdue  University 

School  of  Electrical  Engineerir^ 

Lafayette,  Indiana 

Attn:  Dr.  Julius  Tou 

University  of  Pennsylvania 

Moore  School  of  Electrical  Engineering 

200  South  33rd  Street 

Philadelphia  4,  Pennsylvania 

Attn:  Miss  Anna  Louise  Campion 

Varo  Manufacturing  Company 
2201  Walnut  Street 
Garland,  Texas 
Attn:  Fred  P.  Granger,  Jr. 

CIA 

Room  2447  Y  Bldg. 

Washington,  P.  C. 

Attn:  Mr.  Eugene  Pronko 

Dr.  Saul  Corn,  Director 
Computer  Center 
University  of  Pennsylvania 
Riiladelphia  4,  Pennsylvania 

Applied  Riysics  Laboratory 

Johns  Hopkins  University 

8621  Georgia  Avenue 

Silver  Spring,  Maryland 

Attn:  Supervisor  of  Tochnical  Reports 

Navy  Department 
Washington,  D.  C. 

Attn:  CDR.  J.  C.  Busby,  Code  W3 
Chief,  Bureau  of  Supplies  A  Accounts 

U.  S.  Naval  Avionics  Facility 
Indianapolis  IP,  Indiana 
Attn:  Litrarian,  Code  031.  2 


(1) 
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National  Aeronautics  and  Space  Administration 
Goddard  Space  Flight  Center 
Greenbelt,  Maryland 
Attn:  Chief,  Data  Systems  Div., 

C.  V.  L.  Smith  (l) 

Federal  Aviation  Agency 

Bureau  of  Research  and  Development 

Washington  25,  B.  C. 

Attn:  RD-375/  Mr.  Harry  Hayman  (l) 

Cornell  Aeronautical  Laboratory  Inc. 

P.  0.  Box  235 

Buffalo  21,  New  York 

Attn:  Systems  Requirements  Dept. 

AE  Murray  (l) 

National  Aeronautics  and  Space  Administration 

Theoretical  Division 
S719  Colesville  Road 
Silver  Spring,  Maryland 

Attn:  Mr.  Albert  Arking  (l) 

Stanford  University 
Electronics  Laboratory 
Stanford,  California 

Attn:  Professor  N.  M.  Abramson  (l) 

Lincoln  Laboratory 

Massachusetts  Institute  of  Technology 
Lexington  73,  Massachusetts 

Attn:  Library  (l) 

Commanding  Officer 

U.  S.  Naval  Postgraduate  School 

Monterey,  California 

Attn:  Dr.  Richard  Dorf  (l) 

Waddell  Dynamics,  Inc. 

5770  Soledad  Road 
La  Jolla  ,  California 

Attn:  B.  L.  Waddell,  President  (l) 

Electronics  Research  Laboratory 
University  of  California 
Berkeley  A,  California 

Atln:  Director  (1) 

Institute  for  Defense  Analysis 
Communications  Research  Division 
Von  Neumann  Hall 
Princeton,  New  Jersey 


(1) 


